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Abstract 

The properties of the states of the alternating parity bands in actinides, 

Ba, Ce and Nd isotopes are analyzed within a cluster model. The model 

is based on the assumption that cluster type shapes are produced by the 

collective motion of the nuclear system in the mass asymmetry coordinate. 

The calculated spin dependences of the parity splitting and of the electric 

multipole transition moments are in agreement with the experimental data. 
PACS: 21.60.Ev,21.60.Gx 
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I. INTRODUCTION 



The low-lying negative parity states observed in actinides and in heaviest known Ba, 
Ce, Nd and Sm isotopes are definitely related to refiection-asymmetric shapes There 
are several approaches to treat collective motion leading to refiection-asymmetric deforma- 
tions. One of them is based on the concept of a nuclear mean field which has a static 
octupole deformation or is characterized by large amplitudes of refiection-asymmetric vi- 
brations around the equilibrium shape In this approach the parity splitting is 
explained by octupole deformation. Another approach [|7|,|8|,PJTI1[] is based on the assump- 
tion that the refiection-asymmetric shape is a consequence of alpha-clustering in nuclei 
TT| , p!2|Jl3| . In the algebraic model [[7i|8|,|9|,[l0[| the corresponding wave functions of the ground 



and excited states consist of components without and with dipole bosons (in addition to 
the quadrupole bosons), which are related to mononucleus and alpha-cluster components, 
respectively. The variant of algebraic model including the octupole bosons in addition to the 



dipole bosons has been applied in [0,|1^] to the description of the low-lying negative parity 
states in actinides. In |jl6|,|l^,|18|,|19|] a cluster configuration with a lighter cluster heavier 
than ^He was used in order to describe the properties of the low-lying positive and negative 
parity states. In both models |^J^,p|,p^] and |16|,[1^,|T8|,|T9| the relative distance between the 
centers of mass of clusters at fixed mass asymmetry is the main collective coordinate for the 
description of the alternating parity bands. 

Nuclear cluster effects are mostly pronounced in the light even-even N = Z nuclei with 
alpha-particle as the natural building block. There is a nice relationship between alpha- 



cluster description and deformed shell model ||TT|. It is known from Nilsson-Strutinsky 
type calculations for light nuclei that nuclear configurations corresponding to the minima 
of the potential energy contain particular symmetries which are related to certain cluster 
structures p0|j21|j2^ . By using antisymmetrized molecular dynamics approach p3|j2^ , the 
formation and dissolution of clusters in light nuclei, like ^°Ne and ^^Mg, are described. The 
idea of clusterization applied to heavy nuclei does not contradict the mean field approach. 
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The coexistence of the clustering and of the mean field aspects is a unique feature of nuclear 
many body system. The problem of existence of a cluster structure in a ground state of 
heavy nuclei has attracted much attention, especially, because of the experimentally observed 
cluster decay |^ . The available experimental and theoretical results provide the evidence for 



existence of fission modes created by the clustering of the fissioning nuclei . Indications 



of clusterization of highly deformed nuclei are demonstrated in [P7i|28 . 

The aim of the present paper is a development of the cluster-type model which provides 
not only a qualitative but also a quantitative explanation of the properties of alternating par- 
ity bands. The description of the excitation spectra, ii^A-transition probabilities (A=l,2,3) 



and the angular momentum dependence of the parity splitting [^,0 are the main subjects 
of this paper. Our model is based on the assumption that the reflection-asymmetric shapes 
are produced by the collective motion of the nuclear system in the mass asymmetry coordi- 



nate [0. The values of the odd multipolarity transitional moments (dipole and octupole) 



are strongly correlated with the mass asymmetry deformation of nucleus. In general, the 
value of the quadrupole moment is related to the degree of the quadrupole correlations 
(deformation) in nucleus. However, the collective motion in the mass asymmetry degree of 
freedom simultaneously creates a deformation with even and odd multipolarities. Therefore, 
the calculations of i?A-transition moments are of interest in the proposed model. The single 
particle degrees of freedom are not taken explicitly into consideration since our aim is to 
show that the suggested cluster model gives a good quantitative explanation of the observed 
properties of the low-lying negative parity states. If it is so, this model can serve as a good 
ground for development of an extended model with additional degrees of freedom. 

It should be noted that the first results of calculations of the alternating parity spectra 
for a few actinides within the cluster model are already presented in the Letter |31|]. Besides 
Ra, Th and U isotopes, in the present paper we present the results of calculations of the 
energies of alternating parity bands in 240,242p^^ 144,146,1483^^ i46,i48Qg^ i46,i48n^_ rj.^^ 

electromagnetic transitions are described in this paper with the cluster model for many nuclei 
and the spin dependence of the intrinsic quadrupole transition moment is predicted for ^^^U. 



Simple analytical expressions obtained for the parity splitting and the spectra of alternating 
parity bands are useful for the estimations. The dependence of alpha-clusterization in 
actinides on the angular momentum is shown for the first time. 



II. MODEL 

A. Hamiltonian in mass asymmetry coordinate 

Dinuclear systems consisting of a heavy cluster Ai and a light cluster A2 were first in- 



troduced to explain data on deep inelastic and fusion reactions with heavy ions [32,33,3^ 
The mass asymmetry coordinate 77, defined as r] = {Ai — A2)/{Ai + A2), {\v\ = ^ if ^2=0 and 
Ai = A), which describes a partition of nucleons between the nuclei forming the dinuclear 
system and the distance R between the centers of clusters are used as relevant collective vari- 
ables . The wave function in r] can be thought as a superposition of different cluster-type 
configurations including the mononucleus configuration with |77|=1, which are realized with 
certain probabilities. The relative contribution of each cluster component in the total wave 
function is determined by the collective Hamiltonian described below. Our calculations have 
shown that in the considered cases the dinuclear configuration with an alpha cluster [t] = rja) 
has a potential energy which is close or even smaller than the energy of the mononucleus at 
|?7| = 1 |2^j31[|. Since the energies of configurations with a light cluster heavier than an a- 
particle increase rapidly with decreasing |?7|, we restrict our investigations to configurations 
with light clusters not heavier than Li (77 = rjn), i.e. to cluster configurations near \ri\ = 1. 
The Hamiltonian describing the dynamics in rj has the following form 

where B^j is the effective mass and U [r], I) is the potential. In order to calculate the depen- 
dence of parity splitting on the angular momentum and the electric dipole, quadrupole and 
octupole transition moments we search for solutions of the stationary Schrodinger equation 
describing the dynamics in 77: 



H^niV,I)=En{mn{vJ)- (2) 

The eigenfunctions \l/„ of this Hamiltonian have a well defined parity with respect to the 
refiection t] — > —77. Before we come to the results of Eq. (0), we discuss the calculation of 
the potential U{ri, I), the mass parameter 5^ and the moments of inertia '^{f]) appearing in 
H. 



B. Potential energy 

The potential U{ri,I) in Eq. ([l|) is taken as a dinuclear potential energy for \ri\ < 1 

U{R, V, I) = + B2{v)-B + V{R = R^, V, I)- (3) 

Here, the internuclear distance R = Rm is the touching distance between the clusters and is 
set to be equal to the value corresponding to the minimum of the potential in R for a given 
7]. The quantities Bi and B2 (which are negative) are the experimental binding energies 
of the clusters forming the dinuclear system at a given mass asymmetry r], and B is the 
binding energy of the mononucleus. The quantity V{R,ri,I) in @ is the nucleus- nucleus 
interaction potential. It is given as 

V{R, r], I) = V,oui{R, V) + Vn{R, v) + VrotiR, V, I) (4) 

with the Coulomb potential Vcouh the centrifugal potential Vrot = Tt?'I{I + 1)/(253(?7, R)) and 
the nuclear interaction Vn- In the realization of the cluster model developed in this paper, 
where overlap of clusters is much smaller than in the model of [1^, the choice of the relevant 



cluster configuration follows the minimum of the total potential energy of the system with a 
cluster-cluster interaction taken additionally into consideration. As the result we describe 
the same nuclear properties as in with configurations of clusters having larger mass 
asymmetry and a smaller overlap. 

The potential V{R,t],I) and the moment of inertia Q{ri,R)) are calculated for special 
cluster configurations only, namely for the mononucleus (|?7|=1) and for the two cluster 



configurations with the a- and Li - clusters as hght clusters, respectively. These calculated 
points are used later to interpolate the potential smoothly by a polynomial. The energies 
of the Li-cluster configurations are about 15 MeV larger than the binding energies of the 
mononuclei considered. Therefore, for small excitations only oscillations in r/ are of interest 
which lie in the vicinity of |77| = 1, i.e. only cluster configurations up to Li - clusters need 
to be considered. The potential Vat is obtained with a double folding procedure with the 
ground state nuclear densities of the clusters. Antisymmetrization between the nucleons 
belonging to different clusters is regarded by a density dependence of the nucleon-nucleon 
force which gives a repulsive core in the cluster-cluster interaction potential. Details of the 



calculation of Vat are given in The parameters of the nucleon-nucleon interaction are 



fixed in nuclear structure calculations Other details are presented in |^ . 

Our calculations show that the potential energy has a minimum at \!]\=ria in 
2i8,220,222,224,226j^g^ 222,224,226^1^]^ isotopcs. In ordcr to demonstrate the dependence of 

the potential on the neutron number, we present in Fig. 1 calculated values of U{7]a,I = 
0) = U{ria) of configurations with an a - cluster taking the long chain of Ba isotopes as an 
example. In the neutron deficient isotopes U{ria) is smaller than zero and an a-clusterization 
is more likely. When the neutron number approaches the magic value of 82, the nucleus be- 
comes stiffer with respect to vibrations in t] and U{r]a) is larger than zero. The appearance 
of two neutrons above shell closure is in favor for an a - clusterization. In this case U{ria) 
drops much and again becomes smaller than zero. Further addition of neutrons increases 
the nuclear stiffness with respect to rj vibrations. 

C. Moments of inertia 

The calculation of the moment of inertia Q{r]) = '^{rj, Rm) needed to determine the 
potential energy at / 7^ has been described in For completeness, we repeat in this 



subsection the most important information. As was shown in [^], the highly deformed 
states are well described as cluster systems and their moments of inertia are about 85% of 
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the rigid-body limit. Following this, we assume that the moment of inertia of the cluster 
configurations with a and Li as light clusters can be expressed as 

^{^)=c,{% + % + m,^Rl). (5) 

Here, Q'^, {i = 1,2) are the rigid body moments of inertia for the clusters of the dinuclear 
system, ci=0.85 p8i|3l| for all considered nuclei and mo is the nucleon mass. 

It should be noted that the angular momentum is treated in this paper as the sum of 
the angular momentum of the collective rotation of the heavy cluster and of the orbital 
momentum of the relative motion of the two clusters. Single particle effects, like alignment 
of the single particle angular momentum in the heavy cluster, are presently disregarded. 

For |?7| = 1, the value of the moment of inertia is not known from the data because 
the experimental moment of inertia is a mean value between the moment of inertia of the 
mononucleus (|r7|=l) and the ones of the cluster configurations arising due to the oscillations 
in rj. We assume that 

^{\r,\ = 1) = C2^^{\v\ = 1), (6) 

where is the rigid body moment of inertia of the mononucleus with A nucleons calculated 
with deformation parameters from |^ and C2 is a scaling parameter which is fixed by the 



energy of the first 2+ or other positive parity state, for example 6^. The chosen values of 
C2 vary in the interval 0.1 < C2 < 0.3. So, in our calculations there is a free parameter C2. 
However, this parameter is used to describe the rotational energies averaged over the parity 
and not the parity splitting studied in this paper. 

D. Mass parameter 

The method of the calculation of the inertia coefficient B„ used in this paper is given 



m 



3^. Our calculations show that 5^ is a smooth function of the mass number A. As a 



consequence, we take nearly the same value of B^=20 x lO^mo fm^ for almost all considered 
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actinide nuclei with a variation of 10%. However, for ^^^Th and ^^'^'^^^Ra we varied B„ in the 



range i?^=(10-20)xl0^mo fm^ to obtain the correct value of Eq{I = 0). These variations of 



Bri lead to better results for light Ra isotopes than those in where the obtained values 
of the parity splitting at the beginning of the alternating parity band are smaller than the 
experimental ones. Using a smooth mass dependence of -B^ we get S^=4.5xl0^mo fm^ 
in the Ba, Ce and Nd region. However, better results we obtain for i?^=3xl0^mo fm^. 

For very asymmetric dinuclear systems, we can use simple analytical expressions to 
establish a connection between the relative distance and mass asymmetry coordinates on 



one side and the multipole expansion coefficients P2 and P3 on the other side |28 



7 TT 



Atc 3 



(7) 



Here, Rq is the spherical equivalent radius of the corresponding compound nucleus. One 
finds 



di] 



12 



(l + r^)^/=' + (l-r^) 



1/3 



X 



:i-3,^)+,(i-,-^^+^)-^^^-(^-^)"^^^ 



(8) 



(l+^)l/3 + (l_^)l/3 

In the actinide region for an a -particle configuration, rj ^0.96 and {d/S^/di])'^ ^11.25. Then 
the mass parameters for (3^ and 77- variables are related as 



^ {dPs/dvfBf,^. 



(9) 



If we take the value of Bp^ = 200h^ MeV ^ known from the literature then Br, ~9.3 



xlO^mo fm^. This value is compatible with the one used in our calculations. 



III. INTRINSIC ELECTRIC MULTIPOLE MOMENTS 



Solving the eigenvalue equation (^), we obtain the wave functions of the positive and 
negative parity states for different values of the quantum number I of angular momentum. 
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These wave functions are used then to calculate transition matrix elements of the electric 
multipole operators by integration over rj. The electric multipole operators for a system of 
a dinuclear shape have been calculated ||28| by using the following expression 



Qx, = \/^^ J p'{r)r'Y,,{Q)dr. (10) 

For slightly overlapping clusters when the intercluster distance Rm is about or larger than 
the sum of the radii of clusters (-Ri + i?2); the nuclear charge density can be taken as a 
sum of the cluster charge densities 

pV) = pf(r)+pf(r). (11) 

Using (jTT|) and assuming axial symmetry of the nuclear shape, we obtain the following 
expressions for the intrinsic electric multipole moments 

Qio = 2Dio = 6^(1 - r,')Rm{^ - (12) 
Q20 = e^(l - v')R'miil - ^)|; + (1 + ) + ^20(1) + Q2o(2), (13) 
Q30 = 4^^- v')Rl{{i - vf^^ - (1 + Vf^) + \Rm{{l - rifQ2,{l) - (1 + V?Q2o{2)l (14) 

where the charge quadrupole moments of clusters Q2o{i) ii = 1,2) are calculated with 
respect to their centers of mass. Effective charges for electric dipole and octupole transitions 
are used in our calculations in order to take the coupling of the mass-asymmetry mode to 
the higher-lying giant dipole and octupole excitations ||4^ effectively into account, which 
are not present in the model. 

The charge to mass ratios Zi/Ai and Z2/A2 are functions of rj. For instance, for \ri\ = l 
(mononucleus) this ratio takes the values 0.3-0.4 for the nuclei considered in the paper. 
For the a-particle this ratio is equal to 0.5. The results for the electric dipole moment are 
sensitive to the dependence of Z/ A on rj. In the calculations we parameterize the Zi/Ai 
ratio in the following way. For rja < < 1, the ratio Z2/A2 for the light cluster takes the 
same value as for the mononucleus. For a smaller value of \ri\, we set it equal to 0.5 as for 
the a-cluster. 
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IV. RESULTS OF CALCULATIONS AND DISCUSSION 



A. Calculation procedure 

As was mentioned in Sect. II. B, our consideration can be restricted to cluster configu- 
rations near |?7| = 1. Then it is convenient to substitute the coordinate r] by the following 
variable 

X = rj — 1 if T] > 0, 
X = T] + 1 if rj <0, 

and to use the following smooth parameterization 

[/(x,/) = ^a2fc(/)a;2^ (15) 

k=0 

of the potential U{r],I) from Eq. (^. This formula contains five parameters a2k{I)- If 
a minimum of the potential is located at |x|=a;Q, four parameters are determined by the 
experimental ground state energy, potential energies U (x, /) for and by 

the requirement that the potential has a minimum at |x| = x^. The fifth parameter og is 
necessary to avoid a fall-off of the potential for |x| > xu because of the negative value of 
qq needed to describe correctly U{xLi,I)- We take the minimal necessary positive value of 
as to guarantee an increase of U for |a;| > xu- The ground state energy is obtained by 
solving the Schrodinger equation. Since the ground state wave function is distributed over 
X, the potential energy at x = is not equal to the experimental binding energy of the 
mononucleus. To reach the correct value of the ground state energy Eo{I = 0) = 0, we can 
vary the potential U{x = 0, 1 = 0). In the majority of cases this procedure leads to the 
value of f/(x = 0, / = 0) close to Eq{I = 0). The variation of is also done in the case of 
light Ra isotopes to obtain Eo{I = 0)=0. Besides the barrier height, which determines the 
stiffness of the potential well at x=Xa, the ground state energy Eo{I = 0) of Eq. (0) depends 
also on the frequency of oscillations in x. This frequency is ruled by the value of the inertia 
coefficient -B^. If the minimum is located at x = (|r7|=l), only three parameters Oq, 02 
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and a4 in Eq. ([T5|) are necessary. Potentials with other parameterizations show almost no 
difference in the description of the parity splitting in the majority of considered nuclei. 



B. Parity splitting 

With Eq. (0) we first calculated the parity splitting for the isotopes of Ra, Th, U, Pu, 
Ba, Ce and Nd for different values of the angular momentum /. The results of calculations 
are shown in Tables 1-5. As is seen from the Tables, they agree well with the experimental 



data p],p|,p|,p|,p|,p|,p|,p|,p|,|52| . The largest deviations of the calculated values from 
the experimental ones are found at low spins in some of the considered nuclei. A good 
description of the experimental data, especially of the variation of the parity splitting with 
A at low / and of the value of the critical angular momentum at which the parity splitting 
disappears, means that the dependence of the potential energy on rj and I for the considered 
nuclei is described correctly by the proposed cluster model. The used value of the inertia 
coefficient -B^ is also important. 

Of course, other effects related to degrees of freedom, which are not included in the 
model, like the alignment of the single particle momenta or interaction with other negative 
parity bands with different K quantum number can contribute as well. However, a general 
agreement between the experimental data and the results of calculations shows that the 
simple cluster model used in this paper gives a firm ground for the consideration of the 
alternating parity bands. 

In the considered nuclei the ground state energy level lies near the top of the barrier 
in 7], if exists, and the weight of the a— cluster configuration (Fig. 2) estimated as that 
contribution to the norm of the wave function which is located at \ri\ < rja is about 5 x 10^^ 
for ^^^Ra, which is close to the calculated spectroscopic factor This means that our 

model is in qualitative agreement with the known a— decay widths of the nuclei considered. 

The spectra of those considered nuclei whose potential energy has a minimum at the 
alpha cluster configuration can be well approximated by the following analytical expression 
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E{I) = ^^I[I + 1], if / is even, 



E{I) = ^j^^I[I + l]+6E{I), if / is odd. (16) 



Here, the parity splitting 6E{I) is given as 



,Ein - '''^''r''^ (17) 

l + exp{boJBoI[I + l]) 



with 



Bn = -r 



( 1 



2 \%r]=l) S3(r/ = 



The quantity Bq describes the change of the height of the barrier with spin /. The moment 
of inertia in Eq. (|16[) is given by the expression 

J(/) = w^{I)%r] = 1) + [1 - w^{I)]^{r] = T]^) (18) 

containing a weight function Wm{I) 

which is the probabihty to find the mononucleus component in the wave function of the 
state with spin / of the ground state band. Since Wm{I) decreases with increasing angular 
momentum, J(J) increases with I in agreement with an experimental tendency. The quan- 
tity Wa{I) = 1 — Wm{I) gives the corresponding probability of the a-cluster component. A 
qualitative derivation of the above analytical formulae is given in the Appendix. The con- 
stants ^{\t]\ = 1) = 0.3 X = 1), w^{I = 0)=0.93, bo = tt MeV-^/^ ^^^^ bi=0.2 MeV'^ 
were obtained by fitting the experimental spectra for the nuclei considered (see Fig. 3). 

These formulae clearly demonstrate that there are two important quantities which prede- 
termine a description of the spectra of the alternating parity bands. They are Ei{I'^ = 1~), 
which is determined by the depth of the minimum of the potential at / = and by the value 
of the mass parameter S^, and Bq, which determines the angular momentum dependence of 
Wm{I), i.e. of J(J) and 6E{I). 
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C. i?A-transitions 



With the wave functions obtained, we have calculated the reduced matrix elements of 
the electric multipole moments Q{E1), Q{E2) and Q{E3). The effective charge for El- 
transitions has been taken to be equal to e'l'^ = e(l + x) with an average state-independent 
value of the El polarizability coefficient x=^0.7 This renormalization takes into account 
a coupling of the mass-asymmetry mode to the giant dipole resonance in a dinuclear system. 
In the case of the quadrupole transitions we did not renormalized the charge 62^^ = e 
although an effective charge of 1.35 e describes the data for actinides better as it is seen from 
the results of calculations. For octupole transitions our cluster model Hamiltonian includes 
the octupole mode responsible for description of the shape variation and deformation of 
the nuclear surface. This is the low-frequency collective octupole mode. However, high- 
frequency isovector as well as isoscalar octupole modes are not presented in the model 
Hamiltonian. For example, to simplify the consideration the charge asymmetry coordinate 
is not independent dynamical one but is rigidly related in our model to the mass asymmetry 
coordinate. The octupole transition operator is not exhausted by the term produced by 
the low-frequency octupole degree of freedom and includes also a contribution of the high- 
frequency octupole modes. For this reason for the octupole transitions the effect of the 
coupling of the low-frequency octupole mode to the high frequency mode should be taken 
into account by the octupole effective charge. The estimate of this effective charge is given 
40| . The combined effect gives ^63'°'^ ~ (0.5 + 0.3r2)e [|^. So, we have taken the effective 



m 



charge to be equal to t^^lroton = l-2e for protons and tj^l^utran = 0-8^ neutrons. 

The results of these calculations are listed in Tables 6, 7 and shown in Figs. 4- 
10. The obtained values are in agreement with the known experimental data for Q^^^ 
4|g|g|g|g3g|g|j5gj51|j5|j5|]. Only in 224Ka and ^^^Ba (for / = 7) the calculated 



values of Dxq are larger by factor of four than the experimental Diq. In Th the isotopic 
dependence of the dipole moment is well reproduced. The higher multipole moments are in 
agreement with the calculations of Ref. [p3]. Taking into account the collective character 
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of our model and the absence of the parameters to fit the data, the description of the ex- 
perimental data is rather good. It should be also noted that the experimental data on the 
dipole moment have some uncertainties. 

The angular momentum dependence of the reduced matrix elements of the electric dipole 
operator is presented in Figs. 4 and 7 for ^^^Ra and ^^^Nd, respectively. The calculations 
qualitatively reproduce the angular momentum dependence of the experimental matrix ele- 
ments The same is true for the reduced matrix elements of the electric quadrupole 
and octupole operators (Figs. 5,6 and 8,9). 

Fig. 10 illustrates the angular momentum dependence of the calculated intrinsic tran- 
sition quadrupole moment. It is interesting that the cluster model shows an increase of 
the quadrupole moment with angular momentum in the transitional nucleus ^^^Ra and a 
constant dependence in the well deformed isotope ^^^U. Staggering shown in Fig. 10 for both 
^^^Ra and ^ss-q nudgj jg explained by the higher weight of the a-cluster component in the 
wave functions of odd / states (see Fig. 2). This cluster configuration has larger quadrupole 
and octupole deformations. 

The calculated results for the E3-reduced matrix elements in ^'^^Nd exceed the experi- 
mental data for transitions to the ground band HH] . This can be explained as follows. The 
experimental E3 matrix elements connecting the negative parity states of ^^^Nd to the f3 
band are unexpectedly large, about 70% of the matrix elements within the ground band 



54 1 . This shows a considerable fractionation of the E3 strength among the K=0 bands. In 



our model the (3 degree of freedom is absent and all the E3 strength is concentrated in the 
transitions to the ground band. The nucleus ^^^Nd is transitional in its collective properties 
between spherical and deformed nuclei and the f3 anharmonicity is quite large. This explains 
a fractionation of the E3 strength between the ground and the (3 bands. The summed E3 
strength for ^'^^Nd corresponds to the intrinsic transitional octupole moment of ~ 2000e 
fm^ (instead of ~ 1500e fm^ for the transitions in the ground band) which agrees with 
the calculated value. 
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V. SUMMARY 



We suggest a cluster interpretation of the properties of the alternating parity bands in 
heavy nuclei assuming collective oscillations in mass asymmetry degree of freedom. The 
existing experimental data on the angular momentum dependence of parity splitting and on 
multipole transition moments are quite well reproduced. This supports the idea that cluster 
type states exist in heavy nuclei. The characteristics of the Hamiltonian used in the calcu- 
lations were determined by investigating a completely different phenomenon, namely, heavy 
ion reactions at low energies. Due to this fact, a predictive power of the suggested model 
is quite high. The proposed analytical expression ([16|) for E{I) can be applied to estimate 
the position of the low-lying states which are not yet measured. The calculated stagger- 
ing behavior of alpha-clusterization can be verified by measuring the angular momentum 
dependence of the width of the alpha-decay. 
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APPENDIX A 

Let us assume that the barrier at x=0 separates two minima of the potential (^: the 
minimum with reflection asymmetric deformation (minimum of the potential at a; = Xa) 
and its mirror image. The nonzero penetration through this barrier lowers the energy of the 
levels with even / with respect to the energy of the levels with odd /. With increasing spin / 
the barrier between two minima, which is equal to Ub{I) = Ub{I = 0) + BqI[I + 1], becomes 
higher and the penetration probability goes to zero. According to standard WKB-analysis 
(with higher order corrections) the transmission probability (per tunneling event) for the 
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potential barrier described by the inverted oscillator with frequency hujb{l) for the energy 
hum{l)/'^ above the potential minimum is given by 



^^^^ 1 , fry UtiI)-hulmiI)/2-.' (^^^ 

1 + exp(27r ) 

and hujb{l) ~ \JUb{I) at fixed Xa- Here, ujm{l) is the frequency of the harmonic oscillator 
which approximates the potential U around the a-cluster minimum. Within the semiclassical 
approximation one can neglect the / dependence of uj^ taking uj„i{I) ~ ^^m{I = 0). This 
frequency then essentially determines the rate oJm/'^ at which the wave packet strikes the 
barrier, so that the effective coherent tunneling frequency, i.e. the shift of the negative parity 
states with respect to the positive parity ones, is given by 

5E{I) = ^P{I). (A2) 

TT 

Assuming that {Ub{I = 0) — \TiuJm) is small compared to TiUhiGi) (this case is realized in nuclei 
considered in the paper), we obtain /ico'm/27r =E{I'^ = For those values of / at which 
Uh{I = 0) is much smaller than BqI{I + 1) we obtain Eq. ([TTD with the fitting parameter h^. 
However, we found numerically that Eq. (|17|) works quite well also at low /. 



The weight Wm{I) of the mononucleus component in the wave function of the state with 
spin I can be expressed through the ratio of characteristic times Tm{I) and t;,(/) which a 
system spends in the minima and at the barrier, respectively 

= wnfi (TV ^^3) 

The mononucleus configuration is located at the top of the barrier. We neglect below the / 
dependence of the Tm{I). Denote by rfe(O) the value of Tb(/) at /=0. At very high angular 
momentum when the barrier height is mainly determined by the rotational energy and is 
equal to BqI{I + 1) the value of T{,(/) can be determined from the time-energy uncertainty 
relation 
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To combine the two limits at /=0 and for J ^ 1, we use the following expression 



nil) 



h/nio) + Ub{i) - hum/2 1 + n{o)Boi[i + i]/n- 



Substituting this result into ([X^), we obtain 

n{0)/{n{0) + Tj 

The last expression can be rewritten as 

where w^{Qi) = n{0)/{n{0) + r^) and bi = ^r5(0)r„,/(r„ + n{0)). 
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FIGURES 

FIG. 1. Potential energy (full circles) of the a~cluster configuration U{ria) = U{r]a,I = 0) as a 
function of the neutron number in Ba isotopes. 

FIG. 2. Calculated probability of the a— cluster component in the wave function of the state 
with spin / of the alternating parity band for ^^"^'^^^Ra and 226,232rp]^_ 

FIG. 3. Comparison of experimental (solid points) and calculated with Eqs. (p!^)-([l9|) (solid 
lines) energies of states of the alternating parity bands in ^^'^'^^^Ra and ^^■^Th. The fitting param- 
eters are the same for all nuclei (see text). Experimental data are taken from [41,42|. 



FIG. 4. Angular momentum dependence of the calculated reduced matrix elements of the 
electric dipole operator (solid curve) in ^^^Ra. The experimental data (squares) are taken from 
PI 



FIG. 5. The same as in Fig. 4, but for the quadrupole operator 



FIG. 6. The same as in Fig. 4, but for the octupole operator 

FIG. 7. Angular momentum dependence of the calculated reduced matrix elements of the 
electric dipole operator (solid curve) in ^^^Nd. The experimental data (squares) are taken from 



FIG. 8. The same as in Fig. 7, but for the quadrupole operator 

FIG. 9. The same as in Fig. 7, but for the octupole operator 

FIG. 10. Angular momentum dependence of calculated intrinsic quadrupole transition moments 
in 226Ra and ^^^\J. 
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TABLES 

TABLE L Comparison of experimental {E^xp) and calculated {Ecaic) energies of states of the 
alternating parity bands in 232-222'p]^_ Energies are given in keV. Experimental data are taken 



from 
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TABLE II. Comparison of experimental (Eexp) and calculated {Ecaicj energies of states of the 
alternating parity bands in ■^^'^^^^^Ra and ^'^'^''^^-^Pu . Energies are given in keV. Experimental data 
are taken from [^,42,43]. For ^^'^'^^^Ra, the parameter C2 was adjusted to the 6"^ state. 
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TABLE III. Comparison of experimental {Eexp) and calculated {Ecaicj energies of states of the 
alternating parity bands in 238-232^^ Energies are given in keV. Experimental data are taken from 
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TABLE IV. Comparison of experimental {Eexp) and calculated (Ecaic) energies of states of the 
ground state alternating parity bands in ^^^"^^'^Ba. Energies are given in keV. Experimental data 
are taken from [^,44]. The parameter C2 was adjusted to the 6"^ state. 
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TABLE V. Comparison of experimental (Eexp) and calculated {Ecaicj energies of states of the 
ground state alternating parity bands in ^^^'-^^^Ce and -"^^^'^"^^Nd isotopes. Energies are given in 



keV. Experimental data are taken from |41,45,4f:,47,48,49|. The parameter C2 was adjusted to the 
G"*" state. 
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TABLE VI. Calculated and experimental intrinsic multipole transition moments. The values 
of the dipole moment Diq are given for those values of the nuclear spin / for which there are 
experimental data. These values of / are shown in the second column. The experimental data are 
taken from @,|l],||Ji|j5|,|53|. 
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TABLE VII. Calculated and experimental intrinsic multipole transition moments for Ba, Ce 
and Nd isotopes. The values of the dipole moment Diq are given for those values of the nuclear 
spin / for which there are experimental data. These values of / are shown in the second column. 
The experimental data are taken from 1 44, 45, 46, 47, 4^, 49, 50, 51 1. 







^10 




^r20lu — > z J 








(e fm) 


(e fm) 


(^e im j 


e im ) 


e im j 


(^e im j 




calc. 


exp. 


calc. 


exp. 


calc. 


exp. 


Da, 




n n7i (^^\ 
u.u ( n^iu ) 


zou 


OZl 


izyo 






U.zuy yi=o) 


n ^ A { Q\ 
U.i4(o j 










oa 


U.UOO [I — oj 


U.U0^4 j 


zoo 


ODO 


11/17 






n 1 70 i' T—7\ 


U.UO 1 \0} 










i48Ba 


0.095 (1=7) 




306 




938 




i46Ce 


0.121 (1=7) 


0.11(2) 


313 


305 


1669 




i46Ce 


0.160 (1=11) 


0.20(2) 










148 Ce 


0.152 (1=7) 




364 


436 


1771 




146Nd 


0.071 (1=11) 


0.17(2) 


264 


276 


1627 




148Nd 


0.115(1=1) 


0.24(6) 


370 


400 


2161 


1500 




0.222 (1=8) 


0.24(3) 











29 



3 I 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 

- Ba 

2 - • 




1 - • 

2 - 

L_? 1 I 1 I 1 I 1 I 1 I 1 I 1 I 1 I lJ 

60 64 68 72 76 80 84 88 92 



N 




Spin (I) 



> 

CD 



LU 



2250 - 



1500 - 




2250 - 



1500 - 



750 - 



2 4 6 8 10 12 14 16 18 20 22 



Spin (I) 



1400 



1200 - 



1000 - 



^ 800 



V 



600 - 



400 - 



200 - 




6500 



6000 - 226Ra 
5500 - 



5000 I- 
^ 4500 
^ 4000 

A 

3500 



3000 - 



O 2500 - 



V 2000 - 





I 




I 



A 
m 



3750 - 
3000 - 
2250 - 



W 1500 \- 
V 750 k 
h 


-500 - 
-1000 - 



A 



-1500 - 



m 

O -2000 - 
^ -2500 - 
-3000 - 







t 



6 8 
I 



10 



12 14 



